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ABSTRACT 

We argue that existing models of the ideal pulsar magnetosphere are incorrect because of the 
improper treatment of the singular current layer outside the light cylinder. We simulated the 
axisymmetric pulsar magnetosphere in the Force-Free limit of Strong-Field Electrodynamics. It 
turns out that even in the Force- Free limit, some field lines enter the singular current layer which 
lies beyond the light cylinder in the equatorial plane. As a result: (i) about 10% of the Poynting 
flux is dissipated between 1 and 1.5 light cylinder radii, (ii) there is no singular current layer 
within the light cylinder. 

1. The pulsar magnetosphere problem 

Consider a magnetized ideally conducting ball 
rotating in vacuum. If the ball rotates fast enough, 
is large enough, and magnetized enough, it will 
create charges around itself. This can happen 
even before the Schwinger field is reached, by 
the tree-level QED processes, and also by puUing 
charged particles from the surface. Thus, the spin- 
ning magnetized conducting ball creates a mag- 
netosphere around itself. It is thought that neu- 
tron star magnetospheres are described by this 
model, and the infinitely rich pulsar emission phe- 
nomenology somehow follows. Then, to under- 
stand pulsar emission, one should first understand 
the large-scale structure of the pulsar magneto- 
sphere. 

The ideal pulsar magnetosphere was "drawn" 
by Goldreich and Julian (1969) and calculated by 
Contopoulos, Kazanas and Fendt (1999). These 
calculations were improved by Gruzinov (2005) 
who also gave the pulsar power formula and 
some exact results for the axisymmetric pulsar. 
Spitkovsky (2006) did the non-axisymmetric case. 

All these results, we think, are incorrect. In 
reality, a fair fraction of the magnetic field lines 
enters the singular current layer outside the light 
cylinder. As a result (for axisymmetric pulsar): 

1. Some 10% of the Poynting fiux is damped 
between 1 and 1.5 light cylinder radii. 

2. There is no singular current layer within the 




Fig. 1. — Blue: poloidal current, rBg, above 
the equator. Red: dissipation rate, E • j, be- 
low the equator, outside the star. Black curve: 
Poynting flux through the cylindrical surface, cx 
/ rd2;(i?e_B2 Be), as a function of radius. Sim- 
ulation parameters (in pulsar units, VL = c = 1): 
rs = 0.5, r„ = 3, = 4.5, a = cTs = 333, 
?7 = 1.5 X 10^3 
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light cylinder. 

The error of all existing calculations comes from 
the improper treatment of the singular current 
layer outside the light cylinder. Here we describe 
simulations of the pulsar magnetosphere which re- 
solve all singularities - we don't use any boundary 
conditions, not even at the surface of the star. 

Our numerical simulations arc done in the 
FFE (Force-Free Electrodynamics) limit of the 
SEE (Strong-Field Electrodynamics). We describe 
FFE and SEE in the next section. In §3, we de- 
scribe the simulations. 

2. FFE and SFE 

Both FFE and SFE arc plasma physics mod- 
els which describe the plasma implicitly (Gruzinov 
2008). Namely, one solves the Maxwell equations 

a.F"^ = -f, (1) 

or in the 3-1-1 split, 

dtB = -VxE, atE = VxB-j, (2) 

supplemented by some Ohm's law, which gives j 
in terms of the electromagnetic field only. 
In FFE, the Ohm's law is 

F''^> = 0, Eo = 0. (3) 

Here the scalar Eq is the proper electric field, de- 
fined by 

B^-E^ = B^-E^, BoEo='E-B, Eo>0.{4) 

The physical meaning of the FFE Ohm's law 
is as follows. For any electromagnetic field, at 
any event, there is a one-parameter family of good 
frames, where E is parallel to B. FFE postulates, 
that in any good frame, the electric field vanishes 
and the current flows along the magnetic field. 

In SFE, the Ohm's law is 

BoF''-j, = EoF''''j,, j^r = -a^El (5) 

Here F is the dual tensor, and a is the conductivity 
scalar. 

The physical meaning of the SFE Ohm's law 
is as follows. At each event, the family of good 
frames contains the best frame, where the charge 
density vanishes, and the current aEo flows along 



the common direction of the electric and magnetic 
fields. If £^0 = 0, the charge density p has to move 
at the speed of light, so that j^Ji^ = 0. 

In numerical simulations, one uses the 3-1-1 
split, and the FFE Ohm's law becomes 

._ (B-VxB-E-Vx E)B + (V • E)E x B 

(6) 

The SFE Ohm's law is 

. ^ pE X B + (p^ + j^a^Eiy/\BoB + ^pE) 
•' B^ + Ei 

(7) 

where 

In the limit of high conductivity, SFE reduces to 
FFE (Gruzinov 2008). This might seem strange, 
because SFE postulates that the 4-current is al- 
ways space-like or null-like, j^jfj, < 0, while FFE 
admits time-like currents, j'^j^ = — > 0. 
But it turns out, that SFE handles the time-like 
currents by constantly switching the direction of 
the null-like current, such that the time- averaged 
current becomes time-like. 

FFE can be applied only to initial electromag- 
netic fields of special geometry - with the electric 
field everywhere smaller than and perpendicular 
to the magnetic field. SFE applies to arbitrary 
initial field. 

FFE is ideal, the electromagnetic energy is con- 
served. SFE is semi-ideal, the electromagnetic en- 
ergy is non-increasing, but it remains exactly con- 
stant for all fields with Eq = 0. 

3. Numerical simulations 

In cylindrical coordinates (r, 6, z), assuming ax- 
isymmetric field, we numerically integrate the fol- 
lowing equations 

a^B = -V X E r]AB, (9) 

atE = V X B-j-Fr?AE. (10) 

The small diffusivity rj is added for regularization. 
The equations are solved in a volume 

r < r„, \z\ < Zm- (11) 

The boundary conditions at are 

drB = drE = 0. (12) 
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The boundary conditions at 

B = E = 0. (13) 
The initial conditions are 

B = E = 0. (14) 

The Ohm's law outside the star 

+z^> rl (15) 

is given by eq. ([7]) , with the field invariants regular- 
ization of (Gruzinov 2008). The Ohm's law inside 
the star is the standard relativistic Ohm's law in 
a moving medium 

j = cr,7,(E + vxB)+v(p-cr,7,E-v)+j,. (16) 

Here Ug is the conductivity of the star; v is the 
uniform rotation, and je is the external current - 
both purely toroidal, inside the star: 

ve = rir, jee oc r; (17) 

7s is the Lorentz factor of v. 

The numerical scheme was just the direct prim- 
itive discretization of the PDEs. The only subtlety 
is the time step. Since SFE handles the time-like 
current regions by permanently switching the sign 
of Bq , one needs to reduce the time step (after sat- 
uration of the fields) to get accurate final results. 

We increased the spatial resolution and con- 
ductivities a and CTs , and decreased the diffusivity 
J], until the final magnetosphere outside the star 
showed clear signs of saturation. The final state is 
shown in Fig.l. As one can see, the magnetosphere 
is different from what it was thought to belli 
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^In fact, the dissipative pulsar magnetosphere simulations of 
Gruzinov (2008b) agree with the results of Fig.l. However, 
the numerical scheme of Gruzinov (2008b) was not good 
enough, and at that time we thought that as the conductiv- 
ity increases, the dissipative pulsar magnetosphere would 
approach the standard FFE magnetosphere. We now be- 

lieve that this assumption is wrong. This 2-column preprint was prepared with the AAS 

macros v5.2. 
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